Separate Universes Do Not Constrain Primordial Black Hole Formation 
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Carr and Hawking showed that the proper size of a spherical overdense region surrounded by a flat FRW 
universe cannot be arbitrarily large as otherwise the region would close up on itself and become a separate 
universe. From this result they derived a condition connecting size and density of the overdense region ensuring 
that it is part of our universe. Carr used this condition to obtain an upper bound for the density fluctuation 
amplitude with the property that for smaller amplitudes the formation of a primordial black hole is possible, 
while larger ones indicate a separate universe. In contrast, we find that the appearance of a maximum is not a 
consequence of avoiding separate universes but arises naturally from the geometry of the chosen slicing. Using 
instead of density a volume fluctuation variable reveals that a fluctuation is a separate universe iff this variable 
diverges on superhorizon scales. Hence Carr's and Hawking's condition does not pose a physical constraint 
on density fluctuations. The dynamics of primordial black hole formation with an initial curvature fluctuation 
amplitude larger than the one corresponding to the maximum density fluctuation amplitude was previously not 
considered in detail and so we compare it to the well-known case where the amplitude is smaller by presenting 
embedding and conformal diagrams of both types in dust spacetimes. 
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I. INTRODUCTION 

Einstein's general relativity (GR) permits the fascinating 
possibility that matter and fields interact with geometry in 
such a way that compact regions form from which temporar- 
ily no information can escape. Furthermore, assuming reason- 
able properties of matter and fields, the appearance of such a 
closed horizon (classically) implies the existence of a curva- 
ture singularity in the interior and that information can never 
escape |3 1. These configurations are black holes (BHs). 

In a spherically symmetric spacetime, which we assume 
throughout the paper, a horizon and hence a BH forms as soon 
as the areal radius 7? of a region which determines the measur- 
able surface area A = AnR^, decreases and becomes smaller 
than twice the enclosed mass Achieving this matter con- 
figuration requires strong compression of the material and in 
most circumstances pressure forces will overcome the gravi- 
tational forces preventing the BH formation. 

A dying star is a well-known exception. Depending on its 
mass gravity can overcome the forces exerted by the neutrons 
in the collapsed stellar core. In this way the quantum physics 
of neutrons determines the minimal mass of a stellar BH to be 
larger than the Tolman-Oppenheimer-Volkoff limit, which is a 
few solar masses. 
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Another possibility is that BHs form in the very early uni- 
verse from piimordial density perturbations |[I]|2l|4]|5l. These 
BHs are called primordial black holes (PBHs) and are the sub- 
ject of this paper (we do not consider the formation involving 
the collapse or collision of defects, see for instance |6-9|). 
The formation process is bottom-up like in the ordinary struc- 
ture formation scenario since only perturbations smaller than 
the Hubble radius R < // ' can collapse. The maximum size 
of a PBH forming in the early universe is approximately given 
by the Hubble radius H^^ at the time fhc when the perturbation 
crosses the horizon (see Sec. [II]). As a consequence, PBH sizes 
range in principle between the Planck length and today's Hub- 
ble radius. The actual mass spectrum of PBHs then depends 
on the spectrum of primordial density perturbations and the 
equation of state. Numerical simulations of PBH formation in 
a radiation filled universe fl&HTSl confirmed Carr's estimate 
1 2 1, that the initial fluctuation amplitude of a density perturba- 
tion must be very high in order not to disperse but to form a 
PBH. 

PBHs would open up windows into many interesting 
aspects of GR. For instance, the formation process in a 
radiation-filled universe is a critical phenomenon ifTJl [141 [T6l - 
18 1. The mass spectrum of PBHs strongly depends on the 
probability distribution of metric fluctuations ^ [191 l20l . 
which depends in turn on the mechanism that generates the 
primordial density perturbations 12JL-26J, such that a (non- 
)observation of PBHs constrains these mechanisms. And most 
importantly, while Hawking radiation 1271 of stellar BHs is 
much colder than the CMB and hence undetectable, PBHs 
could have akeady evaporated or evaporate at the cuiTent time 
|l2]|28l[29l, which could give rise to a gravitational wave back- 
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ground f30'-'32'l. Additionally, quantum gravity effects may 
alter the evaporation process at its final stage such that there 
could have remained a Planck mass remnant for each evapo- 
rated PBH ll33l l34l . while PBH evaporation could have cre- 
ated the baryons Il35l437l . Many more reasons for studying 
PBH and more references can be found in ll38l . 

While there is still no direct observation of a BH, there is 
convincing evidence for the existence of stellar BHs |,39J and 
supermassive BHs in the centers of galaxies Il40ll4n . in par- 
ticular our own Galaxy Il42ll43l . In contrast there is no trace of 
PBHs and their existence is strongly constrained by observa- 
tions I44ll45l . However it is still possible that both dark matter 
and supermassive BHs may be the results of PBH evaporation 
1331 and PBH clustering |46 47|. 

The intention of the first part of this paper is to clarify a cer- 
tain aspect of PBH formation which originated from the ear- 
liest treatments of PBHs 1 1 , 2 , 5 1 and can still be found in re- 
cent works on PBH formation, e.g. in |[T5l [TSl l23l |48H50| and 
probably many more. Carr and Hawking considered the col- 
lapse of an initially superhorizon-sized homogeneous spheri- 
cal overdensity surrounded by an otherwise flat FRW universe 
|[T]|2l (see Fig.[T]). They estimated the necessary amplitude of 
the fluctuation for PBH formation using an energy argument: 
the gravitational energy at the moment of maximal expan- 
sion, where the kinetic energy and the expansion rate is zero 
in the interior, has to overcompensate the internal energy den- 
sity that will cause the pressure forces during the subsequent 
contraction. For a radiation-filled universe, they found that 

— 1/2 

this requires the overdensity to have a size Qm 7 where 
ftn - g{tm) is the density inside the fluctuation at maximal 
expansion fm. On the other hand, it was noted that the space 
inside the fluctuation at can only have positive curvature 
determined by This leads to a maximal proper size of 
the overdensity since the region would otherwise close up on 
itself and represent a separate universe (SU). They concluded 

— 1/2 

that having no SU requires ^ Qm 1 1 1- Hence Harada and 
Carr concluded that a fluctuation must be finely tuned in order 
to become a black hole but not a separate universe |49|. One 
aim of this paper is to show that the no-SU condition actually 
does not pose a constraint on density fluctuations, and hence 
fine tuning is not necessary. We see this from several points 
of view in Sees. [ll]- 



IV Carr's and Hawking's no-SU con- 



dition will be quickly rederived in Sec. |ll] keeping numerical 
factors in order to show why it does not constrain fluctuations. 
Defining the amplitude of a fluctuation in terms of a curvature 
fluctuation variable ( in Sec. [Hi] we will conclude that a SU 
corresponds to ^ = oo on superhorizon scales and hence SUs 
are physically impossible and irrelevant for PBH formation. 
In contrast to a result in |2|, we find in Sec.|IV|that the density 
fluctuation 6 features a maximal value 6mnx which is indepen- 
dent of the no-SU condition and instead is a consequence of 
the spatial geometry of the fluctuation determined by the slic- 
ing and the Hamiltonian constraint. Furthermore it will be 
clear that this maximum arises naturally if one expresses 6 
in terms of To every 6 (except 5max) belong two diff'erent 
^'s; one is larger and the other smaller than ^^ax = ^(<Jniax)- 
We will point out the relevance of the rather conceptual find- 
ings concerning SUs for the existing PBH literature in Sec.[V| 




FIG. 1. Sketch of an idealized fluctuation: flat FRW for r > vi, and 
closed FRW for r < r„ 



The second part of the paper is devoted to gaining intuition 
about PBH formation: in Sec. VI we compare the collapse 



of two spherical overdensities with the same 6 but different 
( patched to a flat dust FRW universe using the LTB metric. 
The comparison will be carried out with the help of embed- 
ding diagrams and conformal diagrams, which enables us to 
intuitively understand the differences and similarities of both 
types of fluctuations. 

Throughout the paper we consider highly idealized models 
of density fluctuations. This is because they are analytically 
tractable and at the same time exhibit all the necessary fea- 
tures to clarify and show some new aspects of PBH formation. 



n. BASIC EQUATIONS AND NO-SU CONDITION 

In a spherically symmetric spacetime inhabited by a perfect 
fluid T - {g + p)n®n + pg the metric g can be written in 
comoving coordinates, where n = 1/A^e,, as 



1 +E' 



-dr^ +R^dQ^. 



(1) 



A prime denotes a derivative with respect to the comoving ra- 
dial coordinate r and dQ^ denotes the metric of the 2-sphere. 
Energy density g, pressure p, areal radius R, lapse and en- 
ergy E are in general functions of r and f. Using metric ([T]) 
in the Einstein equations yields the so-called Misner-Sharp 
equations |51 1, which can be conveniently written using the 
Misner-Sharp mass 



M = 4k [ gR^R'dr, 
Jo 



(2) 



that is the gravitating mass enclosed by r. In general, M does 
not coincide with the integral of g over the proper volume 
since 

dVB = AnR^R'dr + 47tR^R'(1 + EY^'^dr = dVj. (3) 

The Misner-Sharp equations are given by the Hamiltonian 
constraint. 



R^ = A^^ 



2M 



(4) 
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the momentum constraint, 

E = - 

the evolution equation. 



l+E p' . 

2 —R, 

Q + p R' 



M = -AnpR^R, 



and the Euler equation, 

(In A?)' 



Q + p 



(5) 



(6) 



(7) 



These equations are rederived in |j521 using the ADM formal- 
ism. In order to close the system of equations, we must pro- 
vide an equation of state f{Q, p) -Q and so as to solve 
we give initial conditions on an initial f=const-slice at fi: 



R(ti,r)-air, g(ti,r) and g(ti,r). 



(8) 



where ai is a constant. This is only one possible choice for a 
set of initial conditions. 

Consider now a spherical overdense region, whose homo- 
geneous interior (r < r^) is patched via a rarefied transition 
region to a flat FRW universe (see Fig.[T]l. We can then choose 
for r > ri, 



N^l, E^O and R^b(f)r. 



(9) 



Let us assume that the size of the fluctuation R\ » i/r' is 
initially much larger than the Hubble radius of the flat FRW, 
where H = bib, such that p' ^ Q and hence from (|5]l and (|7| E 
is approximately time independent and approximately con- 
stant The interior then evolves like a closed FRW universe 
as long as this approximation is valid. This means that for 
r < r„. 



N 



-r and R ^ a{t) r . 



(10) 



Requiring p > -g/3, the horizon // ' will grow faster than 
R{t, r„) and overtake it around the time when the size of the 
overdense region becomes maximal. At this time Rsn vanishes 
and we get from Q 



2M„ 

Rm 



(11) 



Since inside the fluctuation g is spatially constant, Eq. (|2]i 
gives M - 4n/3gR^. Inserting this together with ( fTO] ) into 
( [TT| i, we obtain 



I6ngn 



6_ 



(12) 



The right hand side is exactly the Ricci scalar of a 3- 
sphere with radius a^. The density of the flat FRW back- 
ground g^-a - gm will not deviate very much from the central 




FIG. 2. Comparison of r wAlX radial coordinates. Both circles show 
an embedded 3-sphere (Q and (p directions suppressed). 



density and from the Friedmann equation (|9| and Q it follows 
that ^ R^, such that gradients become important and the 
approximation ([TOjl breaks down at least near the matching re- 
gion when the fluctuation starts to collapse. Let us introduce 
a new, approximately comovin^ radial coordinate ;^f defined 
by 



dr 



Vl +£ 



(13) 



This change of coordinates will overcome any issues at points 
where E - -\ and enables us to treat "normal" PBH forming 
fluctuations f < ^max, the marginal d - Jmax and fluctuations 
with ^ > ^tnax within the same coordinate system. For a closed 
FRW universe, where E - -p-, introducing the ^-coordinate 
solves the problem at r = 1: 



[-^^ + r^d^ = [dx^ + ■i\x?xd^^) , (14) 



which was just a coordinate singularity. The above metric is 
the metric of a 3-sphere with radius a, and r parametrizes this 
3-sphere by projecting down the areal radius r of the equato- 
rial unit ball (ranging from to 1) and so the upper hemi-3- 
sphere is not covered (see Fig.|2]l . The new radial coordinate 
X covers the whole 3-sphere and corresponds to its latitude 
(running from to tt). 

From this we observe that the proper radius L = axa of the 
overdense region is constrained 

L < Lsu = • (15) 
Using the Hamiltonian constraint Eq. (fT2]l gives at fm 



i'SU - 




(16) 



An equality - Lsu would mean that the assumed over- 
density was in fact never part of our universe but instead a 



^ Note that the negligence of gradients is sometimes called "separate uni- 
verse approximation". But this has nothing to do with the notion of a sep- 
arate universe used in this paper 



^ That X is approximately comoving is a consequence of the approximate 
time independence of E. 
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FIG. 3. Embeddings of spatial slices of different types of density fluctuations (angular dimensions suppressed). The left figure shows a case 
where ^ < <fn,ax (, Xa < ^Z'^' S < (5max) Called type I, the middle figure shows the case ^ = <fn,ax (, Xa = 'r/S, 6 = iSmax) and the right one shows a 
case where f > ^^ax LXa > ^/2, S < S^.^) called type II. 



SU. This no-SU condition ensures that there should be no 
matter configuration that is topoloeically impossible. Or in 
the words of Carr and Hawking [ I n 

The - = ^ttT''*^ constraint equation implies 

that the 3-geometry of this hypersurface has positive curva- 
ture of order fi in the region where the rate of expansion is 
zero. If this positive curvature extended over a sufficiently 
large region, the spacelike hypersurface would close up 
on itself to form a disconnected compact 3-space of radius 
about fjT^^^. In this case the region would form a separate 
closed universe which was completely disconnected from 
our Universe. Such a situation would not correspond to a 
black hole formed by collapse of matter in our Universe. This 
shows that for black hole formation pR^ < 1. 

To obtain the estimate in the last line, which relates the 
areal radius R„i instead of the proper radius Lm to we 
could argue that in a positively curved space 

R<L (17) 

at r 9^ and hence ( [T5] l is violated if ( fTSl l is violated. 

< (18) 

\ 8em 

Although our derivation of the no-SU condition ([TSj seems to 
be correct and coincides with the cited result of [1 1, the last 
step ( [T7] i relating proper size L and areal radius is a poor 
estimate. Looking at the embeddings of spatial slices at f„ 
for fixed am and different in Fig. [3] we see that the limit 
Ljn Lsn actually requires R^n — > instead of a violation of 
([TSj. From the spherical geometry (see Fig.[3]l it is also clear 
that 

R < /?max = a (19) 



Their fi is our g. 



is the maximum value R^ can take for fixed Am and thus 

R,^<J-^. (20) 
\ 87r£)m 

It is then obvious that the no-SU condition cannot be vi- 
olated because the spatial geometry ^19\ of the overdensity 
ensures the stronger inequality ( pOj i. It is important to note 
that there is a huge conceptual difference between the situa- 
tions when an observable like is constrained by avoiding a 
topology change ( [18] ) or by simple geometry ( |20] i. 

A further important remark is in order concerning the gen- 
erality of our result. We used idealized fluctuations and ne- 
glected the pressure effects of the transition region because 
our argument is purely geometrical. It does not depend on 
Einstein's equations nor on the specific shape of the overden- 
sity. In order to talk about a "radius" let us assume that the 
overdensity is approximately spherically symmetric, which 
should be a good approximation in case of a high amplitude 
fluctuation L53J . For an inhomogeneously positively curved 
region 1541 . there exists an analog of inequality ( [T5] l, where 
a is replaced by the maximal curvature radius of the region. 
However since ( fTT] ) still strictly holds we see that R is actu- 
ally constrained by the areal radius of the "equator" (analog to 
Eq. ( [T9] l) of the positively curved region and not by the proper 
radius L which is constrained by avoiding a SU. 

Note that we would also obtain the bound ( pO] i if we were 
to take the coordinate singularity at r = 1 in metric ( [T4] i se- 
riously for some reason. At this point E = -I and this was 
believed to indicate a SU IfTSi l55l . This may be the reason 
why type II fluctuations are excluded in the PBH forming pa- 
rameter space in ISOl and why S - 5max should indicate SUs 
|2|. We will see in Sec.[IV]why d^ax corresponds to the up- 
per bound in Eq. ( |20] i and hence does not signal a SU. The 
appearance of a dmax is related to the fact that the gravitating 
mass is given by M ~ Vb instead of M ~ V3 (see ([2| and ([3])). 
The existence of type II fluctuations (x > n/l) is explicitly 
mentioned in ifTTI . but not considered in their PBH formation 
simulations since these fluctuations should collapse anyway 
and are thus not interesting for calculating the PBH threshold 
Xbh- Note also that for an asymptotic Minkowski space these 
configurations are known as semiclosed universes I56J . 
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Another reason why one could beheve that 6 - 5niax marks 
the boundary to SUs is that fluctuations with ( > f^ax actually 
become in a certain sense separate during their evolution on 
sub-Hubble scales. Following the time evolution of such a 



type II fluctuation in Sec. VI we will observe that the neck that 
connects the bubble to the flat FRW gets stretched and tears 
itself apart. One might then be tempted to call the pinched off 
bubble a "separate universe". But as we will see in Sec. VI 



this fascinating process is simply an alternative mechanism 
for PBH formation. But it should be clear from the excerpt of 
|1 1 and it was also pointed out in |49| that a region can never 
evolve into a SU; it must be given as initial condition. 



III. CURVATURE FLUCTUATION 

A physical mechanism that is at least in principle able to 
create type II fluctuations is single field inflation (see (57] 
for an introduction to inflation). To qualitatively see how this 
works, let us first define the curvature fluctuation ( for a new 
radial coordinate s that makes the spatial metric y manifestly 
conformally flat; 



(21) 



and 



The scale factor b is the same as in (|9| if we set r 
^ = for r > ft,. 

Roughly speaking inflation ends when the field which 
drives inflation with its potential V, reaches a certain value 0i 
where both slow roll conditions break down. For the moment 
let us slice spacetime such that is spatially constant which 
makes in the </ir-slice the reheating radiation energy density 
and also the expansion rate approximately constant. Nearly 
all the fluctuations are then in the local number of e-folds 

which are directly related to fluctuations in the value of the 
local volume element 



: e 



3i(i,(,) 2 



sin d-. 



where y is the determinant of the spatial metric y ( |2T| . In this 
slicing the interpretation of how the quantum fluctuations of 
create the metric fluctuations is as follows. Regions where 
the inflaton took more e-folds to reach 0i have a larger vol- 
ume compared to the average. Let us estimate the volume 
fluctuation looking at Fig.|3] It is determined by the ratio of 
the proper volume V-i of the overdense region and the corre- 
sponding average flat space volume Vb one would measure if 
the fluctuation were absent. The proper volume is the volume 
of a 3-sphere cut off at a latitude using ([3]l gives 



V^ix) = 27Tb (x- sin;t'cos;i'), 



(22) 



which corresponds to yfy integrated over the fluctuation. The 
volume of flat space is given by the ball at latitude is: 



47r , , 
VbCt) = sinV ■ 



(23) 
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FIG. 4. Comparison of ^, f and 6 



From this we get an averaged fluctuation ( defined by 

?4ln|^. (24) 

For a marginal fluctuation with Xa = Xmax = we obtain an 
average volume fluctuation ^^^^^ ^ 0.3. We can compare this 
^ to the ^ defined by ( |2T| ). In Appendix |a] we derive that as 
long as the size of the fluctuation exceeds the Hubble radius 
the central volume fluctuation ^(0, f) is approximately given 
by: 



-2 In COS 



(25) 



Fig. [4] shows a comparison between ^ and ^. The fluctuation 
corresponding to ^^ax is ^max - 0.7. 

Both ^ and ^ diverge for;^' ~* ^ ™d hence a fluctuation is a 
SU if and only if on superhorizon scales 



c < 



^su 



(26) 



is violated. Obviously it is impossible to violate this no-SU 
condition. 

It is also important to note that during matter or radiation 
domination, ^ will be approximately constant on super Hubble 
scales as long as spatial gradients are small (see eqs. ( |A6[ ) 
and (|5]l or [58 1) and so the shape of the fluctuation will not 
change much. This implies that a type II fluctuation cannot be 
formed dynamically on super horizonscales during matter or 
radiation domination and must be given as an initial condition 
on superhorizon scales. 

Inflation can produce these initial conditions at least in prin- 
ciple. In the desired comoving region the scalar field <p has to 
make quantum jumps up its potential hill a sufficient number 
of times such that inflation ends in that region on average 0.3 
e-folds later. This should in principle be possible, although it 
may be very unlikely. However, this is a different question. 
All we want to emphasize is that the laws of physics allow the 
creation of arbitrarily large bubbles like the one on the right 
in Fig.[3] 

If inflation took place it created one large bubble connected 
to a mother universe by a comparably tiny neck: the large bub- 
ble is (or hosts) our universe while the neck has the size of the 
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initial false vacuum region in the chaotic or phase transition 
undergoing preinflationary universe. Of course the curvature 
of this bubble does not have to be positive. 



In Sec. VI we will study the time evolution of a type II fluc- 



tuation: Once a neck gets smaller than the horizon of its sur- 
rounding mother universe it will start to stretch and eventually 
pinch off', separating the baby universe from the mother uni- 
verse, leaving behind a black hole in both universes. Why 



and how this happens we will also see in Sec. VI Hence every 
neck attached to a decelerating universe will sooner or later 
be destroyed. The only ambiguous sign left is the PBH. It 
is ambiguous because one can never figure out whether the 
PBH was formed by a type I or II fluctuation (Fig.[3]l since a 
type-I-PBH with a sizexa < ?r/2 has the same mass as a type- 
II-PBH withn-Xn and hence they are indistinguishable for an 
observer in the mother universe. This ambiguity stems from 
the nonmonotonicity of R(x)- The density fluctuation variable 
6 which we will discuss next inherits this ambiguity too. 



IV. DENSITY FLUCTUATION 

In (2], Carr derived from the no-SU condition in the form 
([TSj that there is a maximal density fluctuation dm^x with the 
property that for S < 5max PBH formation is possible, while 
5 ^ <Jmax should indicate a SU. We know already that ([TSj 
cannot be violated because of (j20|. We will now show that 
^raax is a consequence of ( |20] i instead of ( fTSj ) and hence reach- 
ing 5n,ax does not indicate a SU. Starting with the Hamiltonian 
constraint Q inside the fluctuation where g is constant we get 



8n 



-gW- + E . 



(27) 



Like Harrison 1551 . Carr |2j| chooses an initial slice whose 
extrinsic curvature inside and outside the fluctuation is equal 
and constant, such that inside and outside the fluctuation 



Inserting this together with R - a^r and E - 
gives 

6 , 

— + 6Hf - Ibngi . 

Outside the fluctuation, Eq. Q becomes at fi 

Stt. 



(28) 
into ( |27| i 

(29) 



Dividing (|29]l by 6m gives 



(30) 



with 



V "1 J 



6{t) = Ait) - 1 A = -. 

Q 



From the geometry of a 3-sphere (Fig. [3]) we have R\{Xa) ^ ci'i 
which results in the following constraint for 6i: 



5i< 



H. 



-1 \ 



RiiXa) ) 



= 6T 



(31) 



Fluctuations which are initially superhorizon-sized have to 
fulfill (Ji < 1 in this slicing. For fixed Hi the density fluctu- 
ation Si hits this boundary for a size Xa = Xm-dx = 7r/2 because 
RiiXmsLx) - fli- Hence fluctuations with amplitude ^^ax saturate 
the inequality, while fluctuations with x < Xmax and x > Xmax 
fulfifl 6i < 5™", since Riixa) < ^i- To see how the no-SU 
condition Li < Lsu - Oin constrains Si we can divide (|30]l by 
JT^ and use the rather weak condition Ri{ra) < Lsu to get: 



Si < TT^S 



.2 cinax 



(32) 



The no-SU condition leads to a constraint J32| on Si which 
again is ensured by the geometry of space such that the 
no-SU condition does not pose a constraint on Si. 

Because (J™^" depends not only on the initial time but 
also on the initial size Riir^), it is more convenient to eval- 
uate S at the time of horizon crossing fhc, which is defined by 
^hc - ^^hc ■ following, we neglect gradients and pres- 

sure eff'ects in the transition region between the closed and 
flat FRW part. We expect our expression for S to be qualita- 
tively correct, since only around the time fhc this approxima- 
tion starts to break down . Let's see what we find for 

S = Aituc) - 1 

in the approximate synchronous gauge (|9]l and ( [TOj i, when 
we choose as initial condition a simultaneous big ban^ To 
this end we use that the energy densities for closed and flat 
radiation-filled FRW universes are ll57l : 



= 



Stt sin^ rj 



(33) 



where the conformal time rj ranging from < rj < n is related 
to f by 

t = «„,(! - cos 77) . (34) 
Any density fluctuation in the synchronous slicing has then 

at maximal expansion j]^ = n/2. For the marginal fluctuation 

Xa = XmRX = 7t/2 Wifll 

_ _ 1 _1 
Ryu — flm — fm — "i^^m 

we see that it must reenter slightly before maximal expansion: 
the < tm, since R^n < . Because for fixed g^ the marginal 



^ Using the constant extrinsic curvature initial condition would lead to cum- 
bersome expressions unless we would stick to small 5 « 1 . 
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fluctuation has the largest possible areal radius, types I and II 
will reenter earlier and we see already that 6 features a max- 
imum value (Jniax and that it is again determined by the ge- 
ometry and not SUs. The calculation in Appendix [B] shows 
that 

(j ~ -L iiin^xa (8 + fiin^Xa) > (35) 

which is plotted in Fig. [4] We see that the density fluctuation 
amplitude vanishes for a SU with ^ — > oo: 

(5su = (5Cr«=^) = 0. (36) 

The interpretation of ([36| is clear: In the limiting case of a SU 
we have to close the bubble in Fig. [3] We get two disconnected 
spacetimes, a flat and a closed FRW. The maximum value of 
( |35| l is reached for f = ^msa- 

5max=50ro=^/2)-0.6. (37) 

We used in this Sec. a radiation-filled spacetime as an exam- 
ple to show two important properties of 6: the existence of a 
maximum independent of SUs and its vanishing approaching 
a SU. These properties hold in general because they follow 
from the behavior of 7? in a positively curved space and the R 
dependence of M in Eq. (|2|. Since 6 was defined at horizon 
entry, we only need to require that k = pig > -1/3. For a 
general ^ > -1/3 an expression forzlm was calculated l|49l . 

V. CONSEQUENCES 

We briefly point out some consequences for the PBH liter- 
ature regarding SUs. Since the findings are rather conceptual 
they will not change any numerical values found so far but 
they will shed new light on them. 

Musco and Polnarev found Smux - 0.65 for spherical fluctu- 
ation shapes including the rarefied transition region (see Fig. 
11 in jTSll ). They discuss that their 6 is approximately S and 
we see that the result from our simple model agrees quite well. 

But we do not agree with the common physical picture of 
5max, namely that it marks the boundary to SUs ^ [15] [TS] 
|23l |48l - l50l |55l . Rather the appearance of is a conse- 
quence of (0 the spatial geometry which depends on the cho- 
sen slicing and (//) the Einstein field equations, in particular 
the Hamiltonian constraint, which together always ensure the 
no-SU condition given in |r|. The impossibility of its viola- 
tion can be better understood in terms of the curvature fluctu- 
ation ^ Eq. p5| ), where a SU is obtained as the limiting pro- 
cess ^ — > oo during which the Hamiltonian constraint can be 
fulfilled. The utility of the no-SU condition given in [Jj, red- 
erived as relation ^TE\ as well as the no-SU given in 121, red- 
erived as relation P2| i is dubious since they have no geometric 
interpretation compared to relations ( |20] i and pT| ) that enforce 
the no-SU conditions anyway. It seems that we need not worry 
about SUs at all. 

The expression for 6 in Eq. ([35j is nonmonotonous in Xa 
and hence 6 is also nonmonotonous when written as a func- 
tion of ^ using ( |25] l. This means that physically distinct sit- 
uations quantified by the volume fluctuation ^ can have the 



same 6. This result is also valid for inhomogeneous and non- 
spherical density fluctuations. A general fluctuation may be 
considered as type II if the local areal radius R becomes non- 
monotonous. In the context of PBH formation this suggests 
not to use S but instead to take ( as fundamental for specifying 
the non-perturbative amplitude of a superhorizon sized fluctu- 
ation. Shibata and Sasaki did this in 1 12] which automatically 
includes these type II fluctuations and avoids any coordinate 
singularities at dmux which marks the boundary between types 
I and II. 

If one asks for the probability distribution function (PDF) 
of fluctuations it is also quite clear that the PDF for 6 near (J^ax 
can in no way be gaussian since it vanishes for S > 5max- The 
PDF for ^ could in principle be a gaussian near ^max- But the 
initial condition for PBH formation from superhorizon-sized 
fluctuations is ^ > ^bh = 0.7 H- 1.2 [ 12] , which requires a 
nonperturbative method to obtain a rehable PDF, as done for 
instance in f26l. 

During the radiative era the PBH formation threshold Sbh 
is not far from (5,nax: numerical simulations ifTSll found §bh - 
0.5 H- 0.65 for different fluctuation shapes. Depending on 
the PDF of the fluctuation random field ^(x) from which one 
could derive a PDF for shape S and amplitude ^ of fluctu- 
ations, the inclusion of type II fluctuations may be relevant 
concerning the PBH mass power law 

Mbh ^K(S)- m) - SBn(S)y^'^ Mho,- (38) 

for type I fluctuations determined in IfTTl [TSll . where M^or = 
47r/3H^^ is the horizon mass of the unperturbed FRW at hori- 
zon crossing and K, y and 6bh depend on the shape S . The 
power law caused by the critical phenomenon diminishes the 
relevance of the most probable PBH forming fluctuations near 
<5bh for the PBH mass spectrum |48 , 59 1. Even if the ac- 
tual PDF renders type II fluctuations completely irrelevant, it 
is nonetheless interesting to note that the scaling law should 
change near (J^ax to approximately 

MBH-(l+5(0)Mhor. (39) 

The mass is now determined by geometry and has the strange 
property that the PBH mass shrinks with increasing We 
will comment on this at the end of the next section. One can 
find evidence for the behavior ( (39] l in ifTSll . We can extract 
this information from their Fig. 1, considering the last data 
point, which is closest to ^max- There we have approximately 
<5 ~ ^BH = 0.1. Taking the value (Jbh - ^bh = 0.5 from Ca 
(Fig. 10 and a = 0), we get 6 = 0.6. Inserting this value into 
Eq. ((39]l nearly coincides with Mbh = 10"-^Mhoi = l-6Mhoi 
read off from the Fig. 1 in [18]. The authors mention that the 
scaling law starts to change at 5 > 0.01 -i-^bh- And this change 
just seems to signal that the critical phenomenon ceases and 
that the "normal" horizon mass ( |39| l takes over near dmax- 

VI. COLLAPSE 

We now know how to interpret ^^ax correctly. To get a bet- 
ter understanding of what is actually going on during the col- 
lapse of a type II fluctuation, we will compare the collapse of 
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a type II and type I fluctuation with the same 6 but dififerent f . 
Instead of going through calculations we will observe the col- 
lapse using embedding diagrams. Figs. [6]to[T9]and conformal 
diagrams, Fig. 22 For calculational simplicity we take dust 
as the ideal fluid with p - 0. For type I fluctuations note that 
if we choose 6 large enough to escape the scaling law ([38]l, 



the PBH formation process shown here will not qualitatively 
differ from the case where pressure is present. We will com- 
ment on the eff'ect of radiation pressure for type II fluctuations 
in Fig. [14] Setting p - Q turns the Misner-Sharp system into 
that of LTB |60|. The Euler equation (|7]) allows us to choose 
N = I everywhere. The momentum constraint Q and the 
evolution equation (|6| make E and M time independent, such 
that using Eq. (|2]), the only degree of freedom is R. A solution 
R{r, t) of Q with synchronous big bang is given by: 



— (cos 77-1), sin 77 - 77 = 1 E(r) < 

E ,' ' ' M 



9M 



1/3 



,2/3 



£(r) = 



(40) 

where r] is conformal time with < t] < In and t is cosmic 
time with < f < oo. We will not consider the case E > 0. 



A. Initial data 

Because we have chosen R{ti = 0, r) = 0, we cannot give 
Qiir) and pi(r) as initial data, but we can equivalently give E 
and M. In order to treat type II fluctuations without coordi- 
nate problems, we will use the now exactly comoving radial 
coordinate x defined in Eq. ( [T3] l. Writing r(x) = 'I'ix) and 
denoting a derivative w.r.t. by ', we have 

E(x) = f '2 - 1 

such that we can give as initial condition for the simulations 

M(x) and (P(x) 

as depicted in Fig. [5] They are chosen to match the following 
requirements: 

• At every time conditions (/) - (v) given in Appendix [C| 
are fulfilled. 

• The mass, size and formation time of both PBH types 
are equal. This just requires to choose for the sizes of 
the overdense regions to fulfill ;t'n - n - x\- Note that 
Xi\ and xi are time independent for dust. 

• The energy densities inside and outside the black dots 
are homogeneous, while the transit region contains a 
matter free zone (the constant parts of M). Between the 
black dots M and <t> are sufficiently smooth interpolat- 
ing functions. 

• £■ < everywhere 



The first point guarantees a consistent and physically sane cal- 
culation. The second point makes comparison between both 
simulations easier. The third point is chosen to allow an an- 
alytic check of key events during the simulation. That there 
is a matter-free zone has only aesthetic reasons and does not 
influence much the dynamics of geometry. The last point en- 
sures that the three-dimensional space can be embedded in an 

To really understand what the embedding diagrams show, 
let us briefly recapitulate the basic notions of embedding, in- 
trinsic and extrinsic curvature, photon trajectories and hori- 
zons. 



B. Embedding and curvatures 

1. Embedding and intrinsic curvature 

At a fixed time the slice through spacetime is a Riemannian 
3-manifold. Because the space is isotropic we can neglect the 
two angular coordinates § and (p to facilitate the visualiza- 
tion of the geometry: we go to the equatorial plane d- - n/l 
and drop the <p direction. We can easily embed the remain- 
ing one-dimensional line into an which then encodes all 
the intrinsic geometry of the full 3-D slice. Rotating around 
the symmetry axis would restore the <p direction. The graph 
can be found in the following way. Let us denote the auxil- 
iary fourth spatial coordinate by w. The metric of the plane is 
given by 

= dw^ + dR^. 

The induced Riemannian metric on the graph w(R) should co- 
incide with the original radial part of the LTB metric ( |A1[ ) 
because the angular part of both metrics R^dQ^ coincide: 



ind 



dw 
dR 



+ 1 



ori 



This can be rewritten as 



-1- 1 



R'' = 



from which follows 



This gives the w-component of the embedding: 



1*'. 



Note that this only works as long as the argument of the square 
root is non-negative, which is the case if £ < 0. Slices or parts 
of slices where E > Q can only be embedded in Lorentzian 
auxiliary spaces. The thick line in the following plots, the 
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FIG. 5. Mass and energy functions. Full lines correspond to type I and dashed lines to type II 



embedding, that can be thought of as the spatial slice itself is 
the graph 

{R(X,t),w(x,t)) 

parametrized by x- The length of a segment of the graph be- 
tween xi and X2 is identical to the proper length of a radial 
line in the slice between andxi- Also the intrinsic curva- 
ture ^^ti of the original slice is just the curvature of the line (- 
6/(radius of fitted circle)^). But intrinsic curvature is only half 
the truth: 



2. Extrinsic curvature 

The real slice is not embedded in a flat but in a curved 
Lorentzian spacetime. And so there is more geometry hidden 
in the extrinsic curvature K = l/2y. Since the evolution of 
the LTB system is completely determined by the Hamiltonian 
constraint which can be written as li6Tl 

= KijK'j -K^ + 167TQ, (41) 

we see that if we were able to visualize KijK'j - we would 
have visualized spacetime. Fortunately this can be done. With 
X = R'/<1>' we get from 

yj = diag{X-^,R-^, R-^ sin"^ §) 

and the extrinsic curvature 

Kij = diag(XX, RR, RR sin^ i?) 

the sUce (and embedding) intrinsic vector field 

where is the unit vector in ;(f-direction, and the rate of ex- 
pansion of shells; 



which are the amplitudes of and K-if. These functions com- 
pletely determine the extrinsic curvature part of the Hamilto- 
nian constraint: 



C. Photons 

To get a better grasp of time and causality, we will follow 
the way of a few radially moving photons. Their trajectories 
can easily be calculated from the null condition: a photon 
moves on geodesies whose tangent is always null. Hence a 
photon reaches every distance in zero proper time. For a radi- 
ally moving photon this means: 

, lR'\^ . R' 
= -dt^ +l—\dx t/f = ±^4f' 

where + applies for outgoing and - for incoming photons. 
This equation is integrated and the result Xphit) is added to 
the diagram as colored point (R(Xph(t), t), wixpi^it), t)). 



D. Horizons 

A PBH inevitably forms as soon as spacetime develops a 
closed region in which outward traveling photons cannot in- 
crease their areal radius |[3l. The outermost boundary of this 
inner-trapped region is the BH apparent horizon. In a spher- 
ically symmetric spacetime this black hole apparent horizon 
is a centered 2-sphere within a spatial slice on which radially 
outward traveling photons have constant R(xbh)- This means 
that the derivative of R in this outward null direction 

vanishes: 

5v^0rBH(f),f) = o. 

Outgoing photons with slightly larger r have as usual a grow- 
ing R along their radial geodesic, while photons with smaller r 
decrease their R along their radial outward-directed geodesic 
and are trapped. Another apparent horizon which we take 
as the outermost boundary of the outer-trapped region is the 
Hubble horizon on which radially inward traveling photons 
have constant R(xh)- This means that the derivative of R in 
this inward null direction 



(42) 



u = e, e;- ep 
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vanishes: 

Note that the notion of an apparent horizon is only local both 
in space and time: Consider an inward traveling photon at 
Xphih) > Xnih) whose R increases. Since the Hubble radius 
grows this photon can later at t2 be inside ;^ph(r2) < t'Hih) and 
reach the center. We will see this happen in the simulation. 
For a BH horizon, as was mentioned in the Introduction, the 
rules are different: once trapped, always trapped. 
Slightly rearranging equation (|4]) we get: 

R Rl.^ R'^ \ 

M^-iR^-E)^-{R^-^.lj 

= I (dvR ■ duR + 1) . (43) 

From this expression for the mass M we can read off the nice 
result that a necessary condition for PBH formation isR < 2M 
hence that matter must be within its Schwarzschild radius. At 
the Hubble radius this condition is also fulfilled, so we have 
to be a bit cautious. Inner- trapped regions have: 

R<2M and ^ < 0. (44) 

Outer trapped regions have: 

R<2M and ^ > 0. (45) 



E. The embedding diagrams 

From the following plots displaying the collapse of a type 
II (upper panel) and type I fluctuation (lower panel) in a LTB 
spacetime we see the following information: 

• The thick (black) line is the embedded 3d-space of con- 
stant time f, where we keep only the radial direction. 
Rotating around the symmetry axis would restore the 
(^-direction. Intrinsic curvature and proper length 
are directly visible. 



• The thick embedding line is colored orange (grey for a 
b/w print), if the region is outer-trapped, red (dark grey 
for a b/w print) if it is inner-trapped and black other- 
wise. The outer boundary between orange and black 
segments corresponds to the Hubble radius while the 
outer boundary of red and black regions corresponds to 
a BH horizon. 

• Beyond the outer small (black) dot the energy density is 
that of a flat FRW while inside the interior small (black) 
dot the energy density is that of closed FRW. Approx- 
imately the small (black) dots denote the boundary of 
the matter-free region in the transit region between flat 
and closed FRW. 

• To elucidate the causal structure we will follow the tra- 
jectories of three test photons which are symbolized by 
large (colored) dots. 



• Tangential arrows show the vector field Kj^ that mea- 
sures how fast space contracts or expands in X' 
direction. The arrows are tangential and outward- 
directed (green) if space expands and tangential and 
inward-directed (blue) if space contracts in ;i'-direction. 

• Arrows in /^-direction symbolize amplitude and sign 
of Hr, which measures how fast /^-shells contract 
or expand. The arrows are horizontal and outward- 
directed (orange) if the shell expands and horizontal 
and inward-directed (red) if the shell contracts. Note 
that the direction of the arrows only symbolizes this in- 
formation. Since the auxiliary embedding space has no 
physical meaning those arrows do not show a vector 
field. The vectors and would be tangential to 
7?-shells with amplitude 
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Since the Hubble horizon is tiny (there is only a small black region in 

the center) the spatial shapes are nearly unchanged at this early time ^ 

compared to its initial values at fj = 0. Both extrinsic curvatures A* 

and Hft are very large (long arrows) and positive {outward-directed ^G. 7. t = 1.1 

arrows) corresponding to the small Hubble radius. The arrows are top: A third blue photon traveling inward appeared. But it is trapped 

actually 30 times longer and were rescaled to fit the paper. The fol- and hence recedes as long as the embedding at its position is orange. 

lowing pictures however show the correct length. bottom: The Hubble horizon overtakes the orange outgoing photon. 

top: Two outgoing photons, the yellow and orange dots (light grey 

and grey for a b/w print) start in the middle. 

bottom: One outgoing orange photon starts in the middle. 
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FIG. 8. f = 2.1 

top: Horizon crossing occurred and the neck starts to stretch: the 
proper distance of the empty region between the black dots grows. 
The extrinsic curvature K gets smaller, but in the neck it remains 
larger since it has to create the negative intrinsic curvature '•^R. 
bottom: A third green photon traveling inwards appeared. 
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FIG. 9. ? = 3.1 

top: The outer- trapped orange region gets fragmented: vertical re- 
gions naturally remain trapped for a longer time: Moving in these 
regions does not alter R much. 

bottom: The orange outgoing photon now overtakes the Hubble hori- 
zon. 
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FIG. 10. / = 4.1 

rop: Note that type II fluctuations with their local minimal R may 
be seen as wormholes according to the definition C of |62|. To pre- 
vent a static wormhole from becoming a black hole one needs exotic 
matter violating a few energy conditions (see |63 |). In an acceler- 
ating universe a dynamical superhorizon-sized wormhole will never 
reenter the horizon and hence does not require exotic matter to be 
stable. In our decelerating dust Universe we would also need exotic 
matter to stabilize the throat once it entered the horizon. Since we 
do not have exotic matter we see that just below the orange photon 
at the neck, a tiny red future trapped region appeared. This means 
that the mass of the wormhole got larger than half its Schwarzschild 
radius. Once this happened the wormhole ceases to be traversable 
(in the sense that one cannot return) and everything below the outer 
boundary between red and black - the black hole apparent horizon - 
is doomed. Z^*^ in the neck is still positive (green) but Hj^ became neg- 
ative in the neck center which turned the outer-trapped region into 
an inner-trapped one (see (|44j). Note also that the formation of the 
trapped region occurred 

bottom: The orange photon leaves the scene. No BH horizon has 
appeared yet. 
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FIG. 11. / = 5.1 



top: The orange photon nearly rests above BH horizon. The trapped 
contracting red region grows due to the stretching of the wormhole 
driven by the positive A* (green), 
bottom: The blue photon will soon reach the center. 
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The visible region is now completely within the Hubble horizon. _io -5 o 5 10 

top: Also the outer-trapped region in the bubble is disappearing. The ^ 
blue and yellow photon bravely march towards inevitable oblivion. 



FIG. 13. / = 7.1 

top: At maximum expansion /„, the sign of Hr changed in the bubble 
{orange — > red), where a new future trapped horizon has appeared. 
This is the essential horizon for a type I fluctuation. In the neck only 
Hn changed sign and A* remained positive {green) such that the sec- 
ond term in \42\ becomes positive and (visibly) reduces the negative 
in the neck (see \A\)). This also leads to a direction-dependent 
contraction/expansion of space. Since in the lower hemisphere A* 
got negative, space contracts there in every direction. 
bottom: Here the sign of and & changed in the whole closed part 
(which is only part of the lower hemisphere). 
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FIG. 14. / = 8.1 

top: In the interior of the wormhole Hj^ (horizontal arrows) is neg- 
ative {red and inward-directed) everywhere and in the neck center it 
blows up, while ''^'R gets smaller since A* {tangential arrows) be- 
comes also large but stays positive {green and outward-directed). 
Any matter inside the neck would get extremely spaghettified: ra- 
dially stretched by the large positive A* but still volume reduced due 
to the large negative Hr. 

bottom: The blue photon marches upwards. Still no sign of a BH 
horizon. Compare the latitude of the future trapped region in the 
bubble of the type II fluctuation. 



FIG. 15. / = 9.1 

top: The wormhole has been pinched off: a singularity formed that 
crushed the yellow photon. It is also quite interesting that the singu- 
larity formed in a matter free region (between the two black dots). 
Note that the vertical separation of the cusps in the embedding has 
no physical meaning. The inner-trapped region in Fig.|10| appeared 
shortly after horizon entry, suggesting that also in the case of a 
radiation-filled universe the formation of a non-central singularity 
cannot be avoided: pressure forces will try to bring material trough 
the wormhole, which is effective only when the throat is within the 
horizon. Once within it will immediately shrink in the angular di- 
rection and stretch in the radial direction. Both processes slow down 
the pressure driven matter transport. Hence we expect that any type 
II fluctuation with with ^ > > ^max will form a noncentral sin- 
gularity, where ;^„c is a new shape dependent threshold. It would be 
interesting to determine <f„c by a proper numerical simulation. 
bottom: The green photon reaches the center. 
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FIG. 16. t = 10.1 

top: The different trapped regions in the baby universe have merged. 
The blue photon reaches the center 

bottom: Comparing the latitude x of the lower black dots with the 
latitude of the BH horizon in the baby universe in the upper panel we 
see that the BH horizon will form soon. 
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FIG. 17. r= 11.1 

top: The cut-off "elastic band" has now relaxed: In the baby universe 

is everywhere positive. 
bottom: The BH horizon has now formed around the black mass dots 
and sealed the fate of the photons and the matter in the overdensity: 
All the matter and the green photon will be crushed. The blue photon 
is lucky. (Here Mathematica colored some parts greyish red which 
actually should be black. Hence the blue photon is really outside the 
trapped region. The same holds for Fig.|18|l 
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FIG. 19. / = 25.1 

Both photons required a long time (twice the age of the universe 
after PBH formation at t-Bu = 12.57) to finally escape the horizon, 
although they are still inside the matter free vacuole: they must be 
strongly redshifted. After fan both slices are now physically identical 
and indistinguishable: a PBH with Schwarzschild radius R^h = 2.6 
and mass Mbh = ^bh/2 inside an expanding vacuole embedded in 
a flat exact FRW universe. The metric inside the comoving black 
dots is exactly Schwarzschild such that the existence of the FRW is 
undetectable using gravity experiments like Kepler law or peri-PBH- 
shift inside the vacuole. This has been proven by Einstein and Strauss 
(61. 



FIG. 18. / = I2.I 

top: The orange photon still nearly rests above the horizon. The 
baby universe is so close to its big crunch that the blue photon will 
not reach the PBH singularity in the baby universe. 
bottom: Note that the singularity in the type II fluctuation in Fig.[T4] 
was announced by a blowup of K (while got smaller). Here it is 
'-^'R that blows up while K remains small. 
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F. Conformal diagrams 

Despite the great insight into geometrodynamics that em- 
beddings can give us, the causal structure of spherically sym- 
metric spacetimes is best visualized in the form conformal dia- 
grams 1 57 1, where one maps the whole spacetime with a finite 
coordinate system 77, with the property that lightlike radial 
geodesies are straight lines with: 

Xii]) - +7] + const. 

The conformal diagrams for flat and closed dust FRWs are 
depicted in Fig. |20] 

The right diagram shows a flat FRW spacetime. The diago- 
nal boundary is the future null infinity, where all photons 
go (they are reflected at the vertical x - ^ line). All spatial 
slices end for — > 00 in called spacelike infinity, while all 
timelike trajectories terminate for f — » 00 in the future timelike 
infinity i^. The horizontal dashed line is the big bang singu- 
larity. The left diagram showing a closed dust FRW has a 
simpler structure, also the original comoving coordinates t],x 
are already conformal coordinates 77, Jp. Photons are reflected 
at = and x - ^ lines, originate in the lower big bang 
singularity and terminate in the upper big crunch singularity. 

Before we look at the conformal diagrams of our two space- 
times let us first discuss the seeming separation of the baby 
universe seen in the embedding diagrams of the type 11 fluctu- 
ation. One might have wondered whether the formation of the 
baby universe corresponds to a topology change of the spa- 
tial slice which is proved to be impossible (see Li .63 |). Note 
that even a "normal" dust collapse in asymptotic Minkowski 
space which looks quite like a type 1 collapse, can be foliated 
in such an unnatural way to look like a type II (see lISTl box 
23.1). But here we did not change the slicing! So where does 
this change in the behavior of f=const slices come from? 

Let us calculate the spatial shape of ri(t = const,;^') in the 
transition region for both fluctuations. We can easily solve 
(|40ji for 77 which gives 



rj = 2n - 



ri(t,x) = 



\E\ 



3/2 



M 



The main;^' dependence is in the first term. Figure 21 shows it. 
The common feature of type I and II is that rjix) interpolates 
between a constant finite rj in the closed FRW part and in 
the flat part (this ensures that there is no big crunch). We see 
that the f=const-slice of the type II fluctuation has a maximum 
while the type I slice does not. Since the mass function M was 
chosen to be constant in the transition region the whole effect 
comes from \E\^^^. The fact that \E\i\ exhibits a maximum at 
X - 2.2 while \Ei\ does not, explains the strange difference be- 
tween type I and II seen in the embedding diagrams. Note that 
the appearance of a local maximum of \E]i\ - |0n'^ - 1| in the 
transition region, is inevitable and not just an artifact of the 
chosen matching, since <t>ii has to change its sign somewhere 
in the matching region (compare Fig.|5]). This means that in 
synchronous gauge every arbitrary type II fluctuation (which 
may be defined in general by the nonmonotonicity of R{x)) 




t — const 



X = n 

FIG. 20. Conformal diagrams for closed {left) and flat {right) dust 
FRW. 

rj(t = const, x) 




2.0 2.1 



FIG. 21. Spatial shape of 77 in the transition region. Type II is plotted 
dashed, type I full (shifted for better comparison) 



has this appearance of the collapse. Since our model fluctu- 
ation spacetimes consist of parts of the spacetimes depicted 
in Fig. |20] smoothly glued together, the shapes of the resulting 
conformal diagrams for type I and II are obvious. A numeri- 
cal construction of such diagrams from the same data used for 
the embeddings is shown in Fig.|22] Most important is that 
type I and type II have the same conformal structure. From 
the conformal diagrams we see: the f=const-slices in the tran- 
sition region cause all the difference we saw in the embedding 
diagrams. While for type I fluctuations the big crunch is first 
reached in the closed FRW part and later in the matching re- 
gion, for type II fluctuations the big crunch is first reached in 
the matching region and later on in the closed FRW part. And 
there is no apparent topology change if one uses the 77-slicing. 

Although the embeddings show only the case of a spheri- 
cally symmetric dust collapse, it is clear that a slightly aspher- 
ical type II fluctuation in a radiation filled universe will also 
feature a non-central singularity formation in any reasonable 
slicing and provided that ^ > (see Fig.p4|. The only pos- 
sibility to prevent this for ^ < is the time dependence of E, 
such that the neck could widen through the pressure once the 
fluctuation is inside the horizon and turn it into a type I shortly 
before PBH formation. If this does not happen, then the PBH 
mass should be approximately given by Eq. (|39]l. However 
only a numerical simulation can give a definite answer. 
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FIG. 22. Conformal diagrams of type II {top) and I {bottom) spacetime. The coordinates fi,x are chosen to keep f) = rj, x = X the exact 
closed FRW part and to have a ^-simultaneous big bang and big crunch. The full spacelike lines are t = const slices and the ones from f = 0.1 
and to? =12.1 coincide with the slices shown in the embeddings. Close before the big crunch in the exact closed FRW aifj = 2n corresponding 
to f = 12.56, we added an additional line at / = 12.55. Above this in the conformal diagram for the type I fluctuation the slice at f = 13.1 
terminating in the singularity can be seen, showing that the green photon still lives at this time. The magnified corner of the type II diagram 
shows the termination of t > 8.15-slices in the singularity coming from the flat FRW. They reappear from the singularity well visibly on the 
left. This elfect corresponds to the seeming topology change we have seen in the embedding diagrams of the type II fluctuation. For the 
type I fluctuation spatial slices coming from the flat FRW only end in the singularity but do not reappear for smaller x. If we had used the 
^-slicing in the embeddings, the evolution of type II would have looked very similar to the one of type I. Hence the type I and II spacetimes 
are very similar. The two dashed timelike lines correspond to the comoving boundaries that separate the exact closed and flat FRW parts from 
the matching region. They are depicted as black dots in the embeddings. The straight, diagonal {colored) null lines correspond to the large 
{colored) photon dots in the embeddings. Note that in the type I diagram the upper part of the dark {blue) photon line is approximately the BH 
event horizon, while for the type II spacetime it is a photon line starting in the left comer which is approximately the orange photon. 
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VII. CONCLUSION 



we obtain from ( |A1| | and ( |A2| i these two equations: 



Type II > ^max) and mai-ginal (6 - 5max) fluctuations 
sometimes misinterpreted as separate universes do form pri- 
mordial black holes. 

The formation process looks completely different compared 
to the type I class < fmax) in synchronous gauge. This is 
only a gauge artifact of the synchronous and related "natural" 
slicings and can be avoided in the quite unnatural 77-slicing 
used in the conformal diagram. The appearance of dmax is 
independent of the no-SU condition formulated by Carr and 
Hawking ( [T8] l. It is impossible to violate this condition due 
to ( |20l ). The actual limiting value to form a SU obtained by 
closing a type II is 5su — » corresponding to ^ — » 00. The 
nonmonotonicity of 6(() suggests not to use the density fluctu- 
ation variable S but instead the curvature fluctuation variable ( 
if one also wants to consider larger fluctuations, for instance in 
a probability distribution function of fluctuation amplitudes. 

For homogeneous overdensities every type I fluctuation has 
a 5- twin in the type II class related by = n - xn- If there 
is no pressure, both fluctuations form a PBH with exactly the 
same mass although the formation process looks very diff'erent 
for a type II. This is related to the fact that the gravitating 
mass energy that one can probe by using Kepler's law is not 
the integral of the density over the proper volume but over 
the volume of the would-be flat space using the areal radius 
R (see Eq. Q). A diff'erence between type I and II arises if 
there is pressure present. We expect the existence of a shape 
and pressure dependent threshold ^nc with the property that 
for f > ^nc > ^max a nou-ccutral singularity forms, while for 
^nc > ( ^ fmax the type II fluctuation relaxes via pressure 
driven matter loss within the horizon to a centrally collapsing 
fluctuation. To find this threshold it is necessary to perform a 
numerical simulation. 



Appendix A: derivation of the curvature fluctuation 

For this we rewrite the spatial metric 

1 +E 

in terms of the radial approximately comoving;i'-coordinate: 

dr' 



VI + E(r') 



Denoting the inverse by <t>(x) - r, R(0(x), t) = R(x, t) and a 
derivative w.rt. by ' we get: 

''^s^-^^dx'+RHx.t)dQ\ (Al) 
From the definition of x, we obtain the relation: 

£Or) = ^''(;r)-i. 

Writing the metric in terms of the 3d-conformally flat and 
noncomoving radial coordinate s: 



—dx - ce^ds R - ce^s. 

0' ^ 



Taking the quotient gives 

R' 
R0 



-dx — dins. 



(A3) 



(A4) 



Let us assume again that the fluctuation is localized such that 
foTx ^ Xh space is flat: E(x) = or <P'(x) - 1 implying 
'^'(x) - X ~ Xb + 't'iXb)- Flatness also requires a spatially 
constant Choosing ^ = for s > st - s(xh) allows us to set 
r - s hence 't>(x) - s (in the flat space) and we get from \A3\ 
R' - b - c. Inside the fluctuation we integrate \AA\ from;^' to 
Xb which gives s(x, t): 



s(x, t) = 0(xh) exp 



R' 
R0 



idx' 



(A5) 



Inserting this result into the second equation of ( |A3[ l we obtain 

a^^ixJit) 

R f^" R' , 

In the case of a homogeneous overdensity patched to flat space 
at A'i - Xa we have fPix) = sin;^', R - a sinx and R' - a six a'- 



dx' 



sin(yfl) sinx' 
In fact, this converges for;t' — » and we obtain: 

^(0,f) = ^ = -21ncos^. 



(A7) 



(A8) 



Note that xt - Xa (cf- Fig.[T]) is only a good approximation 
as long as the perturbation exceeds the horizon. Otherwise 
gradients in the matching region between Xa and Xb become 
important and hence E and;t'a and;^f^, themselves become time 
dependent. 



Appendix B: derivation of the density fluctuation 

In order to obtain the ratio /Ihc = p/fi at fhc we use the time 
dependence of these two densities determined by the Fried- 
mann Eq. ( |33| ) of a flat and closed radiation-filled FRW, re- 
spectively. Using fm = Om, we compute the ratio at fhc: 



(A2) 



ehc sin'* ;7hc sin"* T/hc 



The condition for horizon crossing is /?hc - or 

dm sin J]hc sinxa = 2fhc sin //he sin^^ = 2v. 

Using ( [34] l to express 77 in terms of f, we obtain 

V2v - sin;^^a = 2v. 
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The solution is 



(if) £ > - 1 to keep the spacetime Lorentzian, hence R' - 
where £ = - 1 . 



4 sin ^Xa + 1 



which we use to obtain 



^(fhc) 



4v' 



4v2 



sin^ ?7hc (2v/ sinxa)"^ 4v2 
1 



= 11 + 4 sin^;^-„ 



and finally 6 becomes 



16 



sinV„ (8 + sin^;^-^) . 



(Bl) 



Appendix C: Conditions on E and M 

In the LTB spacetime the functions M(r) and E(r) can be 
chosen as initial conditions on a given spacelike slice t = const 
in order to model a desired g or spatial geometry. E and M 
cannot be chosen freely; 

(/) E and M must be sufficiently smooth since they enter 
Einstein equations and at r = we need from Eq. Q, 



e.g. E 
r = 0. 



r and M ~ r , since naturally R 



r near 



(///) Q > everywhere, e.g. from Eq. (|5]l M' < whenever 
R' < 0. 

(iv) There is no shell crossing: shells with larger r must not 
overtake shells with smaller r since at those points R' - 
and hence g diverges. 

(v) It is important to note that R' = is not forbidden but 
requires at this point M' - 0. In order to keep the met- 
ric Lorentzian at points where R' = we need there 
additionally E - -\. 

Point (iv) does not have that much physical significance. 
We made an overidealisation: spherical symmetry and 
noninteracting particles. In reality shocks will form when 
the calculation says shell crossing and the solution cannot be 
trusted once a shell crossing has occurred. Note that shell 
crossing even arises for perfect fluids |52|, implying that not 
the missing pressure is the problem but rather the spherical 
symmetry or the perfect fluid description. For a detailed 
treatment of no-shell-crossing conditions see [651, and for a 
complete treatment of points (/ - v) and alternatives see Il66ll . 
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